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ABSTRACT 

We investigate the process of self-regulated star formation via photodissociation of 
hydrogen molecules in low metallicity clouds. We evaluate the influence region's scale 
of a massive star in low metallicity gas clouds whose temperatures are between 10 2 and 
10 4 Kelvin. A single O star can photodissociate H2 in the whole of the host cloud. If 
metallicity is smaller than about 10 -2 ' 5 of the solar metallicity, the depletion of coolant 
of the the host cloud is very serious so that the cloud cannot cool in a free-fall time, 
and subsequent star formation is almost quenched. On the contrary, if metallicity 
is larger than about 10 -L5 of the solar metallicity, star formation regulation via 
photodissociation is not efficient. The typical metallicity when this transition occurs 
is ~ 10 -2 of the solar metallicity. This indicates that stars do not form efficiently 
before the metallicity becomes larger than about 10~ 2 of the solar metallicity and 
we considered that this value becomes the lower limit of the metallicity of luminous 
objects such as galaxies. 

Subject headings: cosmology: theory — early universe — galaxies: formation 
— H II region — ISM: clouds — stars: formation 
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1. Introduction 

After the recombination era, little information is accessible until z ~ 5, after that 
we can observe objects such as galaxies and QSOs. On the other hand, the reionization 
of the intergalactic medium and the presence of heavy elements at high-z suggest that 
there are other populations of luminous objects, which precedes normal galaxies. Thus, 
a theoretical approach to reveal the formation mechanism of such unseen luminous 
objects is very important. 

The formation process of a luminous object is roughly divided into three steps, 
formation of cold clouds by H and/or H2 line cooling, formation of the first generation 
stars in the cold clouds, and the star formation throughout the clouds. However, the 
mass of the first generation stars are estimated through detailed investigation to be 
fairly large flNakamura fc Umcmura 199S , Omukai fc Nishi 1998| ). Thus, the third 



step is disturbed by the feedback from the massive stars formed in the clouds. The 
main feedback consists of two different processes, UV radiation from the stars and 
energy input by supernovae (SNe). Through ionization of H ( Lin fe Murray 1992j ) 



and dissociation of H2 ( [Silk 1977| , |Omukai Sz Nishi 1999j ), ultra violet (UV) radiation 
has negative feedback on further star formation in the host clouds. Especially, H2 is 
dissociated in such a large region that the whole of ordinary low mass cosmological 
object is influenced by one 05 type star ( Omukai Sz Nishi 1999| ). For the case of a 



metal-free gas cloud, the influence region is much wider than the HII region and the 
metal-free host cloud lacks coolant and cannot cool. Thus, next generation stars are 
hardly formed before the first generation stars die. However, the life time of massive 
stars is much shorter than the cosmological time scale and they die as SNe. By these 
SN explosions, the cloud's gas is often dispersed before significant amount of total gas 
is transformed in stars (e.g., Mac Low fc Ferrara 199S| , Ciardi et al. 199£ , Nishi &: Susa 



1999 ). On the other hand, if the gas binding is not disrupted, next generation stars 
are formed in a cloud which is slightly polluted by heavy elements. Even in the case 
that the host cloud is disrupted by SN explosions, if the remnant gas does not escape 
from the host pregalactic object, next generation clouds, which is slightly polluted, 
will be formed and subsequent star formation will follow. In these polluted clouds, 
heavy elements will become important coolants, if their abundances increase to some 
degree. After the host cloud is enough polluted that star formation regulation by UV 
radiation is not efficient, the effective star formation can start. Thus, the pregalactic 
object, which is a cloud complex will evolve into a luminous object such as galaxy. 

In this paper, we investigate the self regulation of star formation via UV radiation, 
and assess the critical metallicity which enables the formation of luminous objects. 
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2. Influence Region of a Massive Star in a Low Metallicity Cloud 

Around an OB star, hydrogen is photoionized, and an HII region is formed. Lin & 
Murray (1992) considered the star formation regulation via photoionization. However, 
the regulation can be efficient outside the HII region via photodissociation of H2 in a 
low metallicity cloud, where H2 line emissions are the most important coolant. 

Although ionizing photons hardly escape from the HII region, photons whose 
radiation energy are below the Lyman limit can get away. Such UV photons 
photodissociate H2, and a photodissociation region (PDR) is formed around the HII 
region. In a PDR in a metal-free cloud, H2 dissociation effect is very efficient so that 
the region which is larger than the whole of cosmological low mass cloud is influenced 
by only one 05 type star ( Pmukai &: Nishi 1999| ). However, after a cloud is polluted by 
heavy elements, the situation becomes complicated, since other thermal processes may 
be important in a PDR in a cloud with heavy elements. In the region, CO molecules 
are dissociated also, since the threshold UV energy of H2 and CO dissociation are close, 
and C, Si and Fe in the gas phase are ionized, since ionizing energies of C, Si and Fe 
are lower than H. Thus, C + , Si + and Fe + cooperate with H2 as main coolants in a low 
metallicity cloud. On the other hand, dust photoelectric heating becomes important 
heating source in a polluted cloud. In this section, we study how much mass in a low 
metallicity cloud is affected by UV photons from an OB star and, as a result, becomes 
unable to cool in a free-fall time. 

To calculate the heating rate (r) and the cooling rate (A) per unit volume, we use 
the rates of Woffire et al. (1995) and the references there in, and the rates of Galli & 
Palla (1999) (processes related to H2). But we do not include the effects of X-ray and 
cosmic ray. We assume the ionization degree x e as 

x e = x H + + x c + + x Si + + x Fe + , (1) 

where Xi is the abundance of the i element. We assume that abundances of heavy 
elements are determined from the cosmic abundances by scaling proportional to z/^qQ 
The adoptive values are xq = W~ 3 - 52 z/z Q , xo = W~ 3 ' u z/z G , xs\ = 10~ 545 z / z Q and 
XFe = 1O _615 z/z0. We assume that all of these elements are ionized. Here we add the 
extra term x H + to evaluate the effect of relic ionization of cosmological recombination 
and/or previous SN, etc.. We investigate for the cases of x H + = 10~ 4 and x H + = 0. 
However, the overall tendency is not affected by the value of x H +- Thus, hereafter, 
we show the results for the case of x-r+ = 1CP 4 mainly. In a PDR, H2 molecules 



The depletion of the gas phase abundance of heavy elements is serious in the interstellar clouds. However, the 
main coolant except H2 is C + and O, and the depletion of C and O is not so large. Moreover, the depletion of C may 
not be serious considering dust formation in SN ejecta (e.g., Kozasa et al. 198£ ) . But it may be possible that the 



depletion of C and O is more significant than the above estimate. In this case, we should consider that the gas phase 
abundances of C and O represent the heavy element abundance, and hence main results in this paper are almost the 
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are dissociated mainly via the two-step photodissociation process by the Lyman and 
Werner (LW) bands photons. For H2 number density, we use the equilibrium value 
with the initial ionization degree ( Omukai fc Nishi 199S| ). This treatment may result 



in overestimation of xn 2 hence in overestimation of cooling rate (see e.g., Nishi fc Susa| 



1999 ). However, to seek the lower bound of the region affected by the photodissociating 
UV radiation from a massive star, we use the equilibrium value. 
H2 is formed mainly via the H~ process 

H + e^ -» H-+7 (2) 
H + H" -» H 2 + e". (3) 

The rate-determining stage of the H~ process is the reaction (^), whose rate coefficient 
k H - is ( |de Jeng 19721) 

k n - = 1.0 x 10~ 18 T s^cm 3 . (4) 
In a PDR, H2 is dissociated mainly via the two-step photodissociation process 

H 2 + 7 -> — 2H, (5) 

where rate coefficient £; 2s tep is given by ( |Kepner et al. 1989| , praine and Betoldi 199^ ) 



fc 2st e P = 1.13 x 10 8 F LW a -1 . (6) 

Here -Flw (ergs s _1 cm _2 Hz _1 ) is the averaged radiation flux in the Lyman and Werner 
(LW) bands. Thus, the equilibrium number density of H 2 under ionization degree x e is 

2 

n H2 = -r x e n u 

™2stcp 

= 0.88 x 10- 26 x e F L -^rn^, (7) 

where kh is the number density of the Hydrogen nuclei. On the other hand, the 
averaged flux in the Lyman and Werner bands is approximately given by 

Aw ^ W 

Then, if self shielding effect can be neglected, nn 2 is proportional to r 2 and only at the 
very outer region, xh 2 becomes abundant enough for efficient cooling. 

However, self shielding effect is important for ordinary clouds considered in this 
paper. If column density of H 2 becomes larger than 10 14 cm -2 , Flw decreases because 
of self shielding (praine and Betoldi 199"6| ). At the outer region where LW band 
radiation is shielded, H 2 is dissociated via thermal collision and rin 2 becomes thermal 
equilibrium value. In this case, H 2 is the more abundant if the temperature is the 
lower. 

Figure 1 shows the change of cooling and heating rates per unit volume (A and T) 
with the distance from the central massive star for the typical cloud (n = 10 cm -3 , 
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T = 3000 K). Here we have assumed the existence of one 05 star, whose mass is 
~ 40 Mq and luminosity of the LW bands is ~ 10 24 ergs s _1 Hz -1 at the center of 
the cloud, and also assumed n, T and x e are constant in space, for simplicity. At the 
inner region where LW band radiation is not shielded, OI and CII line cooling is the 
dominant cooling process, and hence the cooling rate is approximately proportional to 
z. On the other hand, at the outer region where LW band radiation is shielded, H2 
line cooling becomes dominant for metal poor clouds. 
We calculate the cooling time 

t^-maz, (9) 

Aeff 

where n and T are the number density and the temperature of the cloud and 
A e ff (= A — r) is the effective cooling rate. Figure 2 shows the change of the cooling 
time with the distance from the center. 

If metallicity is slightly high (z/zq ^ 10 -2 ), the cooling time is shorter than the 
free-fall time % (= ( 3 2G^m H n ) 1 ^ 2 — ^ x 10 ? y r ( n /10 cm -3 ) -1 / 2 ) at the all region. 
Here G is the gravitational constant, fi is the mean atomic weight and mn is the 
hydrogen mass. Considering HI region, almost all Hydrogen is atomic and fj, ~ 1.4. On 
the contrary, if metallicity is very low, tg < t coo \ at the inner region and tg > t coo \ at 
the outer region. For the low density (n ~ 1cm -3 ) and high metallicity (z/zq ~ 1) 
case, there exists no cool region where net cooling is negative (heating occurs) at the 
inner region. 

For all cases shown in Fig. 2, i coo i < is achieved at the region outer than a 
certain distance from the center, and we call this transition radius as the cooling radius 
(''cool)- However, in the cases of z/zq = 10 _1 and 1, the Stromgren radius, rs, is larger 
than r coo i- If r < rs, the gas is fully ionized and the temperature becomes high, star 
formation is strongly suppressed. Thus, we take as actual r coo i the largest between 
r coo i and rs- In a region r < r coo i, next generation stars are hardly formed before the 
death of the central star, since cooling is inefficient (t coo \ > tg). Thus, we consider this 
region as a influence region. For the very lower metallicity case, z/zq ^ 10~ 2,5 , r coo i 
depends on z very weakly, since main coolant becomes H2 at the outer region. If we 
adopt x H + = 0, bh 2 depends a little on z so that the influence radius depends a little 
on z. But the overall tendency does not change from the case of xh+ = 10~ 4 . 

The Stromgren radius is 

^Mi^r /3 (^M 1/3 - <io) 



2 Note that although the total luminosity of a star depends strongly on the mass, the dependence of the luminosity 
in the LW bands depends is rather weak. 
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where Q* is the flux of ionizing photons by a OB star and Q t ~ 5.1 x 10 49 s 1 for an 
05 star. The mass within r§ is 

M s ^2xlO 4 M ( m X 1 ( rV (11) 

° VlOcm-V \5.1 x 10 49 s-V V ' 

M s is somewhat smaller than the typical cloud Jeans mass (see Figs. 3-5). Note that 
rs and depend strongly on the mass of the central star. Thus, if we consider a less 
massive central star (M ~ IOMq), r s becomes much smaller but other features of these 
figures hardly change. 

For some cases, especially for the high temperature and very low metallicity case, 
tcoo\ can not be shorter than iff for any distance case because of the insufficiency of 
xh 2 - In this case, r coo i is estimated to be oo. 



3. Star Formation Regulation in Low Metallicity Clouds 

To evaluate the strength of the star formation regulation, we calculate 
the ratio of the cooling radius to the Jeans length, r coo f/rj, in the n- 
T plane for the clouds with various metallicity (0 < z/z & < 1). Here, 
rj = ( r* k J n ) 1/2 ~ 2.1 x 10 2 pc(T/3000 K)V 2 ( n /10 cm" 3 )- 1 / 2 . As shown in 
Fig. 3, for the case of z/zq = 10 -2 , n-T plane is divided into two regions. In the 
region with higher density and/or lower temperature, r coo \/rj is smaller than unity. In 
other word, influence radius of a massive star is smaller than the typical scale of the 
cloud (rj). Thus, the regulation is considered to be ineffective. In the other region, 
r cooi/ r J > 1 and the regulation works well. For dense clouds, in the right lightly shaded 
region, i coo i < iff for the whole HI region. In this case, r coo i = rs <C rj. For low density 
and high temperature clouds, in the upper left shaded region, i coo i > iff for the whole 
HI region. In this case, r coo i 3> rj. 

Figure 4 shows r coo i/rj on the n-T plane, for the case of z/z Q = 10 -15 . In this 
case r coo i/rj < 1 for almost the whole region. Thus, star formation regulation by 
UV radiation is inefficient. On the contrary, as shown in Fig. 5, for the case of 
z/zq = 1CP 2 ' 5 , r coo \/rj > 1 for almost the whole region and the regulation works very 
well. The transition occurs when the metallicity is 10~ 2 ' 5 ^ z/z@ 5s 10 -1 ' 5 , and the 
typical metallicity when the transition occurs is estimated as z/z Q ~ 10~ 2 . 

For the extremely low-metallicity case (z/zq ^ 10~ 4 ), the effect of heavy elements 
on the thermal process is almost negligible. Thus, considering the thermal process, it 
is difficult to distinguish a cloud with z/z & & 10~ 4 from a cloud with z/z & = 0. 
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4. Discussion 

For the primordial gas clouds, if we consider the self-shielding effect, the mass 
within the region of influence is obtained as flOmukai fc Nishi 1999 ) 



= 8 x 10 6 Af Q ( -5Lj V V a LLW , f — ) f -^L_ 

u V10- 4 / V10 24 ergs s-iHz- 1 / V3 x 10 3 K / VlOcm" 3 

(12) 

This mass is larger than the typical cloud mass (see Fig. 3-5). On the contrary, for 
the low metallicity gas clouds, the strength of the self-regulation changes with the 
metallicity. 

As shown in the previous section, in the case of very low metallicity clouds 
(z/zq 10~ 2 ' 5 ), star formation regulation by UV radiation is very effective. In this 
case, star formation rate is very low, since only one massive star can stop the evolution 
of the whole host cloud. If SNe do not disrupt the gas binding, the host cloud is 
polluted by heavy elements little by little and the following continuous star formation 
is possible. Even in the case that the host cloud is disrupted by SN explosions, if the 
remnant gas does not escape from the host pregalactic object, next generation clouds, 
which are slightly polluted, will be formed and subsequent star formation will follows. 
After metallicity becomes high enough (z/zq 10~ 2 ), effective star formation can 
start. Thus, the lower limit of metallicity of luminous objects is roughly estimated 
at z/zq ~ 10~ 2 . By the way, a SN release several xMq of heavy elements and the 
typical cloud mass (~ prj) is ~ 10 6 Mq. Thus, if we consider that the cloud mass is 
10 6 Mq and the mass of the heavy elements released by one SN is 4M , the metallicity 
increase per one SN is ~ 4 x 10 -6 and the change of the metallicity is ~ 2 x W~ 4 zq, 
and hence before efficient star formation begins, about 50 cycles of star formation 
and SN explosion is required. This implies that at the effective star formation epoch, 
there should exist some amount of heavy elements and dust and they are scattered 
well. Then, it is expected that it is hard for Ly a photons to escape from the clouds. 
Moreover, the inefficiency of star formation in low metallicity clouds may explain the 
G dwarf problem (e.g., |Rocha-Pinto fc Maciel 1996 ) and may also result in that the 



reionization epoch of the universe may become later than the previous studies (e.g., 



Fukugita & Kawasaki 1994 |Haiman fc Loeb 19971 , |Gnedin fc Ostriker 1997|) . 



The lifetime of very massive stars is ~ 3 x 10 6 yr. However, as noted above, the 
mass dependence of the luminosity in the LW bands is rather weak, and hence even 
if the mass of the central star is smaller (^ 10 Mq), in this case the lifetime of the 
star is longer, star formation regulation can be efficient. Moreover, these less massive 
stars may not evolve into type II SNe. Thus, after the death of the first OB star, 
star formation could occur somewhere in the cloud, and another OB star could form 
successively. Therefore, star formation regulation becomes often serious even for clouds 
whose tg is longer than ~ 3 x 10 6 yr. 

Lin & Murray (1992) considered the regulation only via photoionization. In such a 
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case, the affected region by an OB star becomes the region within a Stromgren sphere. 
For the very lower metallicity case, rg is somewhat smaller than r coo i. However, if 
metallicity is fairly high z/z Q JSJ 10~ 15 , r coo i becomes rg for almost the whole region in 
n-T plane (see lightly shaded region in Fig. 4). Thus, for the higher metallicity case 
(z/zq ^ 1CP 1 ' 5 ), it is good approximation to consider that star formation regulation 
occurs only via photoionization, but in this case influence region is much smaller than 
the host cloud hence the regulation is not effective. 

Recently, it has been shown that earlier phase of the chemical evolution of the 
Galaxy can be explained by the models with the assumption that SN-induced star 
formation is the only star formation process (e.g., [Tsujimoto et al. 199S , Ishimaru 



& Wanajo 1999). As shown above, when the metallicity is lower (z/zq ^ 10 -2 ), star 
formation can occur only after previous massive stars have died. If these massive stars 
dye with SN explosion, we can consider that SN-induced star formation is the only star 
formation process, since stars can form only after SN explosion. 
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Fig. 1. — Change of the cooling and heating rates per unit volume (A and T) with the distance 
from the central massive star for the typical cloud (n = 10 cm -3 , T = 3000 K and z = 10~ 2 z & ). 
The central star is assumed to be one 05 star. PE-h and PE-c in the figure denote the grain 
photoelectric heating and cooling rates, respectively. The others denote radiative cooling rates of 
H2, OI, CII, Sill and Grain, respectively. 
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Fig. 2. — Cooling time i coo i at the various distance from the center for the typical cloud 
(n = 10 cm" 3 and T = 3000 K). The metallicity is z/z Q = 0, 10 -4 , 10~ 3 , 10~ 2 , 10" 1 and 1. 
The dotted line denotes the Stromgren radius, r§ for a 05 star. The dashed line denotes the 
free-fall time, . 
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Fig. 3. — The ratio of the cooling radius to the Jeans length r coo i/rj for the case of z/zq = 10 -2 . 
The numbers in the figure denote the value of r coo i/rj of the contour. For dense cloud, in the 
right lightly shaded region, i coo i < tg for the whole HI region. In this case, r coo i = r§ <C rj. For 
low density and high temperature cloud, in the upper left shaded region, t coo \ > % for the whole 
HI region. In this case, r coo i 3> rj. The values of the Jeans mass are also shown. For reference, 
considering the cosmological objects, related virialized redshift is also shown on the upper axis for 
the flat universe with ^5 = 1. 
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Fig. 4. — The same figure of Fig. 3 but for the case of z/z & = 10 L5 . 



